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Abstract. In this paper we compute the leading order of the Casimir energy for a 
free massless scalar field confined in a sphere in three spatial dimensions, with the 
Dirichlet boundary condition. When one tabulates all of the reported values of the 
Casimir energies for two closed geometries, cubical and spherical, in different space- 
time dimensions and with different boundary conditions, one observes a complicated 
pattern of signs. This pattern shows that the Casimir energy depends crucially on 
the details of the geometry, the number of the spatial dimensions, and the boundary 
conditions. The dependence of the sign of the Casimir energy on the details of the 
geometry, for a fixed spatial dimensions and boundary conditions has been a surprise 
to us and this is our main motivation for doing the calculations presented in this paper. 
Moreover, all of the calculations for spherical geometries include the use of numerical 
methods combined with intricate analytic continuations to handle many different sorts 
of divergences which naturally appear in this category of problems. The presence of 
divergences is always a source of concern about the accuracy of the numerical results. 
Our approach also includes numerical methods, and is based on Boyer's method for 
calculating the electromagnetic Casimir energy in a perfectly conducting sphere. This 
method, however, requires the least amount of analytic continuations. The value that 
we obtain confirms the previously established result. 



1. Introduction 



The Casimir effect is one of the most interesting manifestations of the nontrivial 
structure of the vacuum state in quantum field theory. This effect appears when 
nontrivial boundary conditions or background fields are present (e.g. solitons). In 
1948 Casimir predicted the existence of this effect as an attractive force between 
two infinite parallel uncharged perfectly conducting plates in vacuum [TJ [2] (for a 
general review on the Casimir effect, see Refs. [3j HI [5], [6], [7]). This effect was 
subsequently observed experimentally by Sparnaay in 1958 [8]. Up to now, most 
experiments on the measurement of the Casimir forces have been performed with 
parallel plates [9] or with a sphere in front of a plane [TO] [TTT [12| [13] . Furthermore, 
the configuration of two eccentric cylinders have both experimental and theoretical 
interest [14J. The majority of the investigations related to the Casimir effect concern 
the calculation of this energy or the resulting forces for different fields in different 
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geometries with different boundary conditions imposed P2] . These cases include parallel 
plates P 12], cubes pH OH HI QH M, 121], cylinders [221I231I231I251I2S1I2ZII2E1I22], 
and spherical geometries [IS, [301 EH E2] • This effect has even been investigated in 
connection with the properties of the space-time with extra dimensions [7J [331 El] , and 
nowadays it is known that the Casimir energy depends strongly on the geometry of the 
space-time and on the boundary conditions imposed [35l [36] . In fact, an interesting 
question is the determination of the conditions under which the forces acting on 
the boundaries for closed geometries are attractive or repulsive in arbitrary spatial 
dimensions [37J M, M> M> El S2] • 

The Casimir effect has many applications in different branches of physics. Perhaps 
the first extensive study of the Casimir effect was in Particle Physics and in connection 
with the development of the bag model of hadrons P [30l [43] gH [45l HH1 [47] SH [49]. 
One of the most fascinating and open problems of theoretical physics is the cosmological 
constant [501 EH E21 E21 [51] • This constant is also a candidate for the dark energy [5T ], [53 ] 
[55] . and the Casimir energy has been studied in this connection [5Tl [531 156] . Moreover, 
the presence of the Casimir forces in many different phenomena of condensed matter and 
laser physics have been established both theoretically and experimentally [571 EH EE 60J . 
In particular, the study of the Casimir effect for massless scalar fields is not only of 
theoretical interest but also has direct relevance to physical systems such as Bose- 
Einstein condensates [EH E21 [63] . 

Calculations of the Casimir energy in spherically symmetric configurations have 

attracted the interest of physicists for many years [30l EH [32l SHI S3, SHI SH]- As 

Boyer [61] first showed, the Casimir pressure exerted by the electromagnetic (EM) field 

on the walls of a perfectly conducting spherical vessel is repulsive (see also [65]). The 

reported Casimir energy and pressure of a scalar field confined in a spherical vessel 

with Dirichlet boundary condition is also reported to be positive [661 Ell EHl EH I7D~1 171]. 

Comparison between the reported values of the Casimir energy for closed geometries 

show that these values depend crucially on three major factors: first the details of the 

geometries (e.g. spheres or cubes), second the number of space dimensions, and third the 

boundary conditions. To be very concrete we have collected all of the reported results 

for various cases and display them in Table (0Q). We should mention that all of the 

results obtained for cubical geometry are exact, while the ones for spherical geometries 

are obtained approximately by various numerical methods. As one can see from the 

Table (Op), for the spherical cases in even spatial dimensions, there always remains an 

unresolved divergent factor. For the case of massless scalar fields with the Neumann 

boundary condition, the Casimir energy is always negative regardless of the details of 

the geometry and the number of space dimensions. Therefore, one might conclude 

that there are no surprises there. However, for the case of massless scalar fields with 

Dirichlet boundary conditions there is a sign factor (— 1) D for the case of cubes [72J, 
and (_i)Floor(D/ 2 -i) for the 

case of spheres [69] , where D denotes the number of space 
dimensions. For the case of EM field inside a perfectly conductor there is a sign factor 
(— 1) D+1 for the case of cube, and no obvious sign factor for the case of sphere. 
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Scalar 


M 


-0.33 


-0.260872-0.044716/s 


4 


EM 


Conductor 


-0.044 


-0.197834-0.033768/s 


5 


Scalar 


V 


-0.0025 


-0.000288 


5 


Scalar 


M 


-0.37 


-0.270281 


5 


EM 


Conductor 


+0.021 


-0.006362 



Table 1. The Casimir energies for the EM field (with boundary condition 
consistent with a perfectly conducting cavity denoted by "Conductor") and 
massless scalar fields (with Dirichlet and Neumann boundary conditions 
denoted by "P" and "jV", respectively). The values are displayed for the 
space dimensions D = {2,3,4,5}, for cubical and spherical geometries. The 
quantity s is a regularization factor which should be taken to zero. The values 
are displayed for a cube with unit volume and a sphere with radius a = 1, all 
in units He = 1. 

Now we concentrate on the case of three spatial dimensions. At a first glance, the 
reported results for the case of a massless scalar field confined in a spherical geometry 
with Dirichlet boundary condition seem anomalous because it is the only case for which 
the two geometries do not have the same sign of the Casimir energy, and this has caused 
a controversy in the literature. As far as this controversy is concerned, L.A. Manzoni 
and W.F. Wreszinski [73], try to justify these results by showing that the Casimir forces 
in both cases are repulsive. It has been claimed that for the EM case, the deformation 
of a spherical shell of radius a into a cubical shell of length L with L « 2a, should 
not change the sign and approximately the magnitude of the Casimir energy, when 
the boundary conditions are unchanged [2TL 1741 175] . We believe that their claim is 
reasonable. However, it is inconsistent with the results for the EM field in five spatial 
dimensions and the massless scalar field with Dirichlet boundary conditions in three and 
four spatial dimensions. Therefore, it is worth studying this problem in more details, 
and here we concentrate on the three dimensional case. 

Although, the Casimir energy for the case of the cube is exactly solvable [74"] . 
the presence of divergences inherent in these sorts of problems require regularization 
and analytic continuation methods [76], and this has been a source of criticism for this 
calculation [TTJ [78] . Since the Casimir energy for the case of the sphere is not exactly 
solvable, it requires the utilization of numerical methods. The use of numerical methods 
for problems which include divergences could be a source of even greater concerns 
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about the accuracy of the final results. When the problem is plagued with multitude 
of divergences, the numerical methods should include delicate and carefully planed 
regularization and possibly analytical continuation techniques. The problem addressed 
in this paper is in this category, and we calculate this case using an alternative method 
which requires the least amount of regularization and analytic continuation schemes. 

There are three general methods for calculating the Casimir energy for the spherical 
case, which as mentioned earlier all include numerical parts. The first method is the 
Green function formalism [32]. In this method one encounters an infinite sum of integrals 
over the modified Bessel functions. It can be shown that each of the integrals converges 
and the sum can be calculated numerically by using the asymptotic behavior of the 
modified Bessel functions for all real D. The plot of the Casimir pressure as a function 
of continuously variable D show that the result is infinite for even dimensions [32] . The 
second method is the zeta function technique. For example, G. Cognola et al [69J use the 
zeta function technique to obtain the Casimir energies for spherical symmetric cavities, 
for different boundary conditions, and fields in D space dimensions. The third method is 
the direct mode summation using contour integration [70] . This method is based on the 
direct summation of the frequencies and the main tool employed is the Cauchy theorem. 
This employment lets the authors sum the zeros of the Bessel functions more easily. 

The method we use in this paper is based on the Boyer's method which was used 
for the EM field in a sphere. An important part of the Boyer's method is to confine 
the system inside a similarly shaped but larger shell, and then to subtract the energies 
of two configurations with the same size outer shells and different size inner shells. In 
this subtraction scheme, most of the infinities automatically cancel without any need 
to use analytic continuations. As a matter of fact in this method, the Casimir energy 
is simply the work done on the system in deforming the initial configuration to the 
one under consideration. Therefore, the quantity just defined has an obvious physical 
interpretation. One can then let the radii of the outer shells and the second inner 
shell go to infinity. We shall henceforth refer to this subtraction scheme as the Boyer 
Subtraction Scheme (BSS). We recently used this method to directly compute the lowest 
order Casimir energy for the EM field inside a rectangular waveguide Ref. [79], and the 
lowest order radiative corrections for the Casimir energy for a scalar field for the parallel 
plate problem in various space-time dimensions [80l El] • 

In this paper we use BSS for a massless scalar field confined in a sphere with the 
Dirichlet boundary condition. The analytic part of our calculation is analogous to the 
one used by Boyer for the TE modes, and is done in section 2. In section 3, we present 
the final part of our calculations which, similar to the Boyer's paper, is numerical. 
However, our numerical calculations differ from that paper. We introduce a method 
to extrapolate the sum of all the zero point energies for a given value of the angular 
momentum £. In order to proceed with the summation over £, we first introduce a second 
extrapolation method to obtain the optimal form of the summand. We then sum over 
all £ by using the zeta function analytic continuation technique to remove the infinities 
and extract the finite part. Our final result agrees with the established value obtained 
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in [SSI EZl ESI EH EE] • In section 4, we summarize and discuss our results. 

2. Analytic Part of the Calculation of Casimir Energy in a Sphere 

In this section we setup all the physical and mathematical machinery for the calculation 
of the Casimir energy for a real massless scalar field with Dirichlet boundary condition 
for a spherical geometry. The mathematical part introduced in this section parallels 
closely the one due to Boyer [64] for the computation of the Casimir energy for the EM 
field inside a perfectly conducting sphere. The similarity between our calculation and 
that of Boyer stems from the fact that the expressions for the Casimir energy for the 
TE mode of the EM field is similar to the massless scalar field with Dirichlet boundary 
condition [6] . The important difference is that the term with I = is allowed in the 
latter case. Moreover, there is no analogue of the TM mode in our problem due to 
choice of boundary condition. As we shall see, this forces us to encounter some new 
divergences. A close examination of Boyer's work show that when both modes are 
present these divergences cancel. The final expression obtained in this section cannot 
be solved analytically, and we solve it numerically in the next section. 

The Casimir energy is defined as the difference between the sum of the zero point 
energies of all the modes in the bounded region considered, and the free case. In order 
to find all of the modes inside a sphere we have to solve the Klein-Gordon equation 
for a massless scalar field, with the appropriate boundary condition. Due to spherical 
symmetry, each energy level Eg will be (2t + 1) fold degenerate. Therefore we have the 
following expression for the total zero point energy, 



where uj£ iS and k^ s are the normal modes frequencies and wave- vectors, respectively, and 
the index s refers to the root number for a given value of t. Since there are an infinite 
number of normal modes of increasingly high frequency, this energy is infinite, as usual. 
However, as mentioned before -Ecas. is obtained by the subtraction of this zero-point 
energy from the analogous one for the free case. An essential part of the Boyer's method 
is to enclose this sphere inside a larger one. Then a similar configuration is considered 
but with a different radius for the inner sphere (see Fig. (CQ)). Then the zero point 
energies of these two configurations are subtracted from each other, and many of the 
infinities cancel each other out. Then the radii of all the spheres except the original 
inner sphere (the one with radius 'a') is taken to infinity. Therefore our expression for 
the Casimir energy is given in Eq. (j2j), 




(1) 




(2) 



Using Eqs. ffT|2l) we have, 




(3) 
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Configuration A Configuration B 

Figure 1. Left figure shows two concentric spheres with radii a and R 
(configuration 'A'). The right figure shows a similar configuration with different 
radius for the inner sphere (configuration 'B'). To calculate the Casimir energy, 
the zero point energies of these two spherical configurations are to be subtracted 
according to the Eq. ([2]). The label T refers to the inner sphere and 'II' the 
outer annular region. In the final step, R — ► oo, while n is kept fixed. 



X 



• oo oo 

y*]k e ,a (a) 9 (Afc/,. 0)) + y~] K ttS (a, R) g [\K e ,s (a, R)) - {a -> R/i]} , 

L s=l S=l 

where k^^ s {a) and K^s (oj R) are the wave- vectors for the inner sphere and annular 
region, respectively, and g(\k( )S (a)) = e~~ Xke ' s ^ is a convergence factor which eventually 
goes to one as A — > 0. Moreover the index S is defined analogously to s, but for the 
annular region. 



The imposition of the boundary condition is explained in Appendix A, where the 
relationships between the root indices, s and S, and the wave-vectors k^ s (a) and 
K^s (a, R), are obtained. Using the results obtained in |Appendix A| and the Euler- 
Maclaurin Summation Formula (EMSF), we calculate the limit R — > oo of terms in 
Eq. (|3]) coming from the annular regions in Appendix B There we show how in this 
limit, most of the infinities coming from these regions cancel due to BSS, and all of the 
remaining finite terms approach zero. The final result is that the two indicated sums 
turn into integrals with no remainder. As explained below in more details, the two 
integrals eventually cancel the infinities coming from their respective inner spheres, and 
their residual infinities cancel each other. We finally obtain, 

he 



^Cas. = Jim lim — J2( 2i + X ) 2 k e , s (a)g(\h, s (a)) 



1=0 



+ / dS K i;S (a, R)g(XK e>s (a, R)) - {a -> R/rj} 



5=1 



, (4) 



where the continuous variable S is associated with the discrete root number Sg(a, K) and 
appears when we use EMSF. The above equation is similar to the expression obtained 
by Boyer for the TE mode, with the only difference being the inclusion of i = in 
this case. It is obvious from their definitions that K^s(a) and Se(a,K) are in one to 
one correspondence. Therefore we can change the variable of integration from S to 
the continuous version of the wave-vectors, which we shall denote by K for brevity of 
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notation, and using Eq. (jA.80 we have, 



E Cas . = lim lim ^ J2( 2£ + !) E M<0*(AA*») " / K g(XK) 



+ 



,IK dSe{R : K) Kg(XK) - {a ^ i?M 



'Jf=Jf/,i(o,fl) 

Rewriting Eq. (jSJ and rearranging terms we have, 



.(5) 



^Cas. = Km \im^f^(2£ + l)\\f^k e Aa)g(XhAa)) - f dK ^fi- Kg(XK) 

R^ooA^O 2 ^ JK=Kt A (a,R) dK 



K=Ke A (a,R) 



dK 



dK 
dse(R, K) 



dK 



Kg(XK) .(6) 



-{a - ^//7} 

' K=K ttl (R/ V ,R) 

As shown in Appendix B Eqs. ( TB.2~|B.3j) . in the limit i? — > oo, the wave- vectors 
K^i(a, R) and K^R/r), R) not only decrease exponentially in but also go to zero as 
1/R. So, the difference between the upper and lower limits of the last integral approaches 
zero and since the integrand is not singular, there is no contribution coming from this 
term to i?Cas .■ Also, the lower limits of the remaining integrals can be extended to 
K = 0. 

We can use the simplest form of the EMSF to make the first term in Eq. ([6]) more 
amenable to computation, 



s=l 



s=l 

OO 



dsk^ s (a) g(Xk^ s (a)) + ^[ke,i{a)g(Xk e>1 (a))] 



f°° ( 1\ d 

+ J ds{^s- [s] - -J—[k e>s (a)g(Xk e , s (a))\, 



(7) 



where [se(a, k)} denotes the Floor Function. The one to one correspondence of se(a, k) 
and ki(a,s) is analogous to the previous case which was for the annular regions. 
Therefore, we can again change the variable of the integration from s to k, and by 
adding and subtracting appropriate terms we can extend all of the lower limits of the 
integrals to zero. We obtain: 



y^ y k e>s (a)g(Xke, s (a)) 



s=l 



°° dk dsei 2 k) kg(Xk)- 
k=o dk 



+-[ke, 1 (a)g(Xk eA (a))] + 



dk 

dk[s-[s] - \)^r[kg{Xk)} 



k=0 



k=k t ^(a) 



k=0 



k=0 
1\ d 



2J dk [ 



dk(s-[s]--)-[kg(Xk)]. 



The last term of Eq. ([8]) can be simplified by noting that the Floor Function [se (a, , 
in the indicated domain, and integration by parts finally yields, 



(8) 
= 



k=ki i (a) 
fc=0 



dk[si(a,k) - [s e (a, k)} - ^ -^[kg(Xk)} 



k=ki tl (a) 
k=0 



dk (st(a, k) - \)-JL[ k 9(^k)\ 



^h,i(a)g(Xk £tl (a)) 



" hAa) dk^S-kg(Xk). (9) 
fc=o dk 
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Finally Eq. flS} 

oo 

ke, s {a)g(\ke 



can be simplified to 

/ % \ f°° IdsAa, k) 
(a)) / " 



oo 

dk 

k=0 



dk 



k <j( A//) ! ( s £ (a, fc) - [s e (a, k)] ~ ^) Jr [kg(^k)} 



.(10) 



Upon substituting the expression displayed in Eq. (ITU1) into Eq. (EJ), the divergent terms 
cancel each other (the first term of Eq. (11 01) and the second term of Eq. (jSJ)). Analogous 
cancelation occurs for the other regions. All of the remaining terms are finite. Moreover, 
we can make the integrals dimensionless by making appropriate changes of variables, 
such as x = ak. We finally obtain 



he r 00 / 1 f°° 

^Cas. = Jim lim - £(2£ + 1) - / dx 



sAx) 



[SAX 



\)^[^g((^/a)x)} 



{a - i?M 



:id 



This expression is analogous to the one for the Casimir energy for the TE mode obtained 
by Boyer for the same geometry, except that £ = is included in our expression. This 
integral cannot be done analytically and we compute it numerically in the next section. 



3. Numerical Evaluation of the Casimir Energy 

Analytic calculation of Eq. (llip seems to be impossible. Therefore, we resort to a 
numerical method. We compute the integrals in Eq. ffTTl) for each value of £ separately, 
and then sum the results as indicated in that equation. However, the numerical method 
used by Boyer cannot be employed, since in our problem we do not have the luxury 
of cancelation of infinities between the TE and TM modes and the constancy of their 
sum. As is apparent from Eq. (llip . the integrand for each £ has an infinite number 
of discontinuities due to the presence of the Floor Function in that expression. After 
determining the precise position of the jumps, the integrations are done separately for 
all parts and then all of the results are summed. The integration is over the continuous 
version of the wave number, which extends to infinity. In order to accomplish this 
numerically we compute the integral up to a cutoff M. For fixed £, the quantity M can 
equivalently be thought of as a cutoff over the root numbers s. Then we compute this 
integral for a series of values of M. The results of the numerical integration, indicated 
by I(£, M), are plotted as a function of M. Then by fitting a polynomial function to this 
plot, the asymptotic value can be easily obtained, and gives us the infinite M limit. On 
the other hand we have to take the limit A — > as indicated in Eq. (|lip . An attempt to 
optimize the accuracy of our results as a function of the relationship between M and A 
has revealed that the best choice is A = 1/M. In Fig. (j2J) we display the computational 
technique just described for the cases £ = {0, 1}. The case £ = also serves a secondary 
purpose, which is a crucial part for the rest of our numerical analysis. For this case the 
frequencies are proportional to integer multiples of ir and the original expression for the 
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5 10 15 20 25 5 10 15 20 25 



0.00 IM/ff 0.00 1M/^ 

Figure 2. In the left figure, we plot the difference between our numerical results for 
I = (the data points and their best fit shown using solid line), and the theoretical 
value Eq &s = (in units Jp = 1) as a function of the cutoff M, as explained in the 
text. Note that, our numerical values asymptotically approach the theoretical value. 
In the right part of the figure we show our numerical results for the case I = 1 (the 
data points and their best fit shown using solid line) . For both figures the asymptotic 
values are obtained by fitting a sum of polynomials (r n ) with n = {—5, 0} and are 
shown using dashed lines. 



Casimir energy, Eq. (OQ), can be computed directly as follows: 

n nc v-^ , «c v-^ sir here t , . hen / — 1 \ . . 

s^E* = yEv. = yE 7 = ¥ a-i)- ¥ y. (i2) 

s=l s=l s=l 

Note that in the last step indicated by arrow we have used the analytic continuation 
of the zeta function [7J [82]. In the left part of Fig. (T5]) we have plotted the difference 
between our numerical results for various values of the cutoff M, for the case £ = 0, 
and the theoretical value ^ (in units lp = 1). As it is apparent from this figure, the 
asymptotic values of our numerical results exactly matches the theoretical one. This 
clearly shows the consistency between our numerical methods and the theoretical results 
obtained by analytical continuation of the zeta function. In the right part of Fig. (T5]), we 
plot the actual numerical values of the Casimir energy for £ — 1, since theoretical values 
do not exist in this case. In Table ([2]) all the results obtained for various values of £, up 
to £ = 21, are shown for various values of the cutoff. Note that the values indicated in 
the last column are the asymptotic values which correspond to the cases in which the 
cutoffs have been sent to infinity. 

It can be shown that the integral ( TTTT) is an odd function of v. An appropriate 
ansatz for the integral, which can be obtained by fitting its asymptotic values reported 
in Table©, is 

/0) = Av + Bv- 1 + CV~ 3 + Dv- b + 0(z/T 7 , where v = £ + \ (13) 

where fiv) indicates the absolute value of the asymptotic values of our numerical 
integration (/(£, M)) recorded in the last column of Table©, and A,B,C and D are 
the unknown coefficients to be determined from the data. To find these coefficients we 
fit an increasing sequence of data points to the above ansatz, and then find the overall 
asymptotic values of these coefficients. To be more concrete, we fit the values of the 
following set of sets of data points: £= {{1, 2, 3, 4, 5}, {1, 2, 3, 4, 5, 6}, {1, 2, 3, 21}}, 
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14 


-7.21742 


-7.24387 


-7.24720 


-7.24831 


-7.24887 


-7.25054 


15 


-7.71268 


-7.74288 


-7.74669 


-7.74796 


-7.74859 


-7.75050 


16 


-8.20764 


-8.24183 


-8.24615 


-8.24759 


-8.24831 


-8.25047 


17 


-8.70229 


-8.74073 


-8.74558 


-8.74720 


-8.74801 


-8.75045 


18 


-9.19663 


-9.23956 


-9.24498 


-9.24680 


-9.24770 


-9.25042 


19 


-9.69067 


-9.73834 


-9.74436 


-9.74637 


-9.74738 


-9.75040 


20 


-10.18441 


-10.23705 


-10.24370 


-10.24593 


-10.24704 


-10.25038 


21 


-10.67784 


-10.73570 


-10.74302 


-10.74546 


-10.74669 


-10.75036 



Table 2. The results of numerical integrations are listed for various values of £ and 
the cutoff M. All of these values are denoted by I(£,M). For the last column, we 
define f(u) = \I(£, oo)|, where v = £ + \ 



and for each set (which contains iV data points) we find the corresponding coefficients 
A N , B N ,C N and D N . 

In Fig. we plot An, B n ,Cm and D N as a function of N . These four graphs clearly 
show that there are asymptotic values for the coefficients which we denote by A^, B^, 

and -Dqo. We had 17 sets of data points for each graph, we fit a polynomial with 
terms N~i with j G {—5, —4, 0}, and obtained the asymptotic values, 

= 0.50000 , = 0.00781 , = -0.00105 , = 0.00034, (14) 

so we have, 

f(v) = A oo v + B O0 v- 1 + C O0 v- 3 + D O0 v- 5 + 0(v)- 7 i > 0. (15) 

Before proceeding with our calculation, it is interesting to mention that f(v) has the 
following reported asymptotic expansion (Eq. (3.2) in Ref. [7Dj), 

f{") = ~ + — ^+ 565 , + 0(u)- 7 £>0. (16) 

JK ' 2 128f 32768z/ 3 1048576z/ 5 y 1 y 1 

Comparing these two expressions, it is obvious that they are equivalent, to within the 

numerical accuracy of our calculation. Now we can substitute our optimal functional 

form obtained for the integral, Eq. (Tl5|) . into Eq. (fTTj) to obtain, 

£ 00 ft 00 

^Cas. = t|5>+^(E WW) ~ {« - ^} 

8=1 £=1 
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Figure 3. The plots of An, Bn,Cn and Djy as a function of N, as explained in 
the text. These four graphs clearly show that there are asymptotic values for the 
coefficients which we denote by A^, -Boo, Coo and Dqo- 



he 

2a 



C(-1)tt - 2A 00 C(-2, 3/2) - 2B oo C(0, 3/2) - 2CU(2, 3/2) 



-2^(4,3/2)) 



{a - R/ V } 



he 

' 2a 
-2L> r 



— 7T 

12~ 



-2A. 



25 00 ( - r 



2C n 



, 7T 



4) 



,7T 



6 



16) 



{a -> -R/77}, 



where in the last step we have used the analytic continuations of the generalized zeta 
functions [H [H2] along with the following shift formula for the zeta functions, 



C(n,m) 



00 

(k + m) 

k=0 



m~ n + (k' + m+ iy 

k'=0 



m~ n + C(n,m+ 1). (18) 



Obviously, in this case the last term in Eq. (|T7|) goes to zero as R — » 00. We have 
checked that we can obtain the same finite result, without resorting to any analytic 
continuation, using the Abel-Plana summation formula. We have separated the £ = 
term from the rest of the summation over £, primarily because the value of the integral 
for I = dose not follow as closely as we like the pattern of other data points. We have 
checked that our shift formula allows us to single out this term, even when we need to 
analytically continue the zeta functions. 

Our final result for the Casimir energy for a massless scalar field with Dirichlet 
boundary condition in a three dimensional sphere is, 



^Cas. = +0.00562 ( - 



(19) 



This result has been obtained numerically and is the same as the previously reported 
result [?, EH EE], which were obtained by different numerical methods, to within 0.2%. 
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4. Conclusion 



In this paper, the Casimir energy for a massless scalar field in a sphere with radius a 
with Dirichlet boundary condition is computed in three spatial dimensions. This energy 
has been computed using the BSS which is based on confining the original sphere inside 
a concentric larger sphere and computing the difference between the vacuum energies of 
two configurations which differ only by the size of the inner spheres. Finally, the radii of 
all the spheres except the original one go to infinity. Our result shows that the Casimir 
energy for a massless scalar field in a cube and a sphere have indeed opposite signs, 
contrary to the case of the EM field. As mentioned in the Introduction, and explicitly 
shown in Table ([1]) , the sign of the Casimir energy for closed surfaces depends crucially 
on the type of field considered, the boundary conditions imposed, and the details of the 
geometry and not merely the topology. To be more specific, the claim that the sign 
of the Casimir energy should not have a significant dependency on the details of the 
shapes of the vessels [2TJ EU ES] is, much to our surprise, not true. 



Appendix A. The Imposition of the Boundary Conditions 

In this Appendix, we show how the imposition of Dirichlet boundary condition on the 
inner and outer spheres relates the root indices, s and S, and the wave-vectors ki >s (a) 
and K^s { a , R), which appear in the main regulated expression for the Casimir energy, 
Eq. ([3]). For the configuration A, we obtain the following condition for the inner sphere 

j £ (ah, s (a)) = 0. (A.l) 

That is ki t3 (a) is the sth zero of jt{ak^ s {a)). The ranges of the values of angular 
momentum and inner root indices are I = {0, 1, 2, 3, ...} and s = {1, 2, 3, ...}. Similarly, 
the following relationship can be obtained by combining the results of the imposition of 
the Dirichlet boundary conditions on the two surfaces confining the annular region, 

n e (aK e!S (a,R)) n t (RK e>s (a, R)) ' 1 ' ' 

where K^s{fl,R) is the Sth zero of Eq. (1A.2I) . The ranges of the values of angular 
momentum and annular root indices are I = {0,1,2,3,...} and S = {1,2,3,...}. In 
order to find the relationship between the root indices and the wave-vectors, we use the 
following relationships between the spherical Bessel functions, trigonometric functions, 
and polynomials in 

xje(x) = sin [x — Ai(x) + cos (x — \t^) Bg(x), 

(A.3) 

xrii{x) = — cos [x — A e (x) + sin (x — ^£it) B^(x), 



where 



A >W ~ Zj (2r)!(< - 2r)! ( > 
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BAx) 



i_ *yr (-i) r (£ + 2r + i)\ 

2x ^ (2r + l)\(t- 2r - 1)! ^ % ' 



(A.4) 



r=0 



Using Eqs. flA.3|A.4l) and Eq. (lA.lj) . we obtain the following nonlinear relation between 
integer s and k^ s (a), 



ake, (a) 



-In = 7rs — arctan 



Be {ah,s (a)) 
A e (ak e , s (a))_ 



which can be simplify to, 

7is = — arctan 



je (ak e>a (a)) 



(A.5) 



(A.6) 



Hi (ake >s (a))_ 

Similarly, we can repeat analogous steps to obtain an equation for integer S as a function 
of K^s{.Q>, R) f° r the annular region, 



tiS = arctan 



je (aK^s (a,R)) 



— arctan 



je {RK,, S (a,R)) 
n e (RK e:S (a, R)) _ ' 



(A.7) 



n e (aK it s(a, R))_ 

Comparing Eq. (|A.6[) with Eq. (|A.T[) we obtain the following relationship between the 
root indices, 

S e (a, R, K) = s e (R, K) - s £ (a, K) . (A.8) 

Note that, in this equation the values of the root indices s, S, and the wave-vector K 
can be thought of as having been analytically continued to non-integer values. 



Appendix B. Implementing the R — > oo Limit 



In this Appendix, we show that in the limit R — > oo, the contribution to the Casimir 
energy, reflected by the last two terms in Eq. ([3]), can be greatly simplified. This is 
due to the fact that, as we shall show, most of the infinities coming from the annular 
regions cancel due to BSS, and all of the remaining finite terms approach zero. Using 
the EMSF, the contribution of the annular region in Eq. ([3]) is, 



J2 [Ke,s (a, R) g (XK e , s (a, R)) - K e>s (R/rj, R) g (XK iiS (R/rj, R))] 

dS [K t (a, R, S) g (XK e (a, R, S)) - K e (R/ V , R, S) g (XK E (R/ V , R, S))] 



5=1 



5=1 



+- [K £>1 (a, R) g (XK e>1 (a, R)) - K iA (R/ V , R) 9 (\Kt,i (R/v, R))] 



r=N 



B 



2r 



d 



2r-l 



^ (2r)! [dS 2 ^ 1 

dSB 2N (S-[S]) 
d 2N 



K e (a, R, S)g(XK £ (a, R, S)) - K e (R/ V , R, S)g(XK e {R/rj, R, S)) 



5=1 



dS 2N 



K e (a, R, S) g (XK e (a, R, S)) - K e (R/ V , R, S) g (XK e (R/rj, R, S)) 
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Here we show that the last three terms in the right hand side (rhs) of Eq. (IB.ip give no 
contribution to the Casimir energy. In fact, in the limit R —>■ oo, all terms in the rhs of 
Eq. ( IB. II) . except for the first integral, cancel each other due to BSS or vanish entirely. 
In order to show the vanishing of these terms, we have to prove two things. First, each 
of these terms decreases rapidly as a function of £, so that when it is summed over £ 
with a pre-factor of 11 + 1 it will not give a divergent contribution. Second, each term 
approaches zero in the limit R — ► oo. 

We start with the second term. For the small values of £, the first root goes to zero 
as R~ l . For large values of £, the first root K~e t i(R/r), R) behaves like, 



K til (R/rj,R) 



where 



v 



1/3 



e ~ 



-2u{/3 u -tanh /3„) 



1.23?/ 



and 



(cosh fly 



(B.2) 



rj-^l + Oiv- 2 ' 3 )]. (B.3) 



and v 
1 



+ 1/2. Therefore, the second term becomes, 



R J 



k t ,i(R) 



e 

R 



a(m^) + £) 2 + . 



(B.4) 



R 



k lil (R) + -) +\(k i)1 {R) + -) + 



R 



0(e/R)-0(e'/R), 



For obtaining the second line we have expanded the exponentials since A eventually 
approaches zero. As shown in Eq. (1B.3I) the numerators e and e' decrease exponentially 
in £, and the entire term goes to zero in the limit R — ■> oo. 

An analogous argument can be used for the third term in the rhs of Eq. (IB. II) . In 
the neighborhood of S — 1, we have the following relationship between the wave- vectors, 

K e (a, R, S) = k e (R, S) + 5 e (a, R, S). (B.5) 

where 5g(a,R,S) decreases exponentially with increasing £ and vanishes as R — > oo. 
Analogous relationship holds for any derivatives of the wave- vectors. In fact, all terms 
which depend only upon the outer radius R (and not on the inner ones a or R/rj) are 
canceled due to BSS, and the remaining terms go to zero as R — » oo. 

Now, we proceed to the fourth term. Here, we show that the remaining integral 
decreases rapidly with increasing £. We rewrite the integral as 

d 2N 



dSB 2N (S-[S}) 



5=1 



dS 2N 



< 



2N 



(2N) 

\B-2N 



(2N) 

\B2N 



(2N) 



rS=S* 


dS 


■ Js=i 




1 rS=oo 


dS 


' Js=s* 




1 rS=oo 


dS 


! Js=s* 





d: 



2N r 



dS 2N 
d 2N 



K t (a, R, S) g (XK e (a, R, S)) - {a -> R/ V } 

K e (a, R, S) g (XK e (a, R, S)) - {a -> R/n} 
[K e {a,R,S) g(XK e (a,R, S))] 
[K e (R/r),R,S) g(\K £ (R/r),R,S))} 



dS 2N 
d 2N 



dS 2N 



(B.6) 
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where S* = S^R/r], R, u/b) and b is constrained by: R/rj < b < R. For large values of £, 
we can repeat analogous arguments given above to show that for the first integral the 
terms of the integrand which depend only upon the outer radius R, are canceled due 
to BSS, and the result of the integration of the remaining terms decrease exponentially 
with increasing I and vanish for R — > oo. As for the second term in the rhs of Eq. (IB. 61) . 
the high order derivatives of the wave-vector (a, R, S) for large value of I are, 

dmKl dS* ,S) =0(^ 1 " m ). for u/R<K t {a,R,S)«:v/a 
= 0(v~ m/3 ), for K e {a,R,S) ~ {u + ru 1/3 )/a 
= 0(u 1 - m ) } for u/a<£K e (a,R,S) (B.7) 

where v = £+1/2 and m > 1. Thus if a sufficient number of terms are taken in the EMSF 
(Eq. (IB. II) ) the high order derivatives of the wave- vector, Kg, decrease with increasing 
v. Moreover, when derivatives operate on the damping factors, giXKt (a, R, S j) = 
e -XK e (a,R,s) ^ eac h derivative extracts a factor of A, 



d 



dS^ 



1 



—g(\K e (a,R,S))~\0 -g(XK e (a,R,S)) . (B.8 



R 



Now, if N in Eq. (IB. II) is chosen to be large enough, each term arising from 
jji^w [Ki (a, R, S) g (XKg (a, R, S))) , decreases rapidly with increasing v or has many 
actors of A. The terms which have factors of A approach zero when A goes to the zero, 
and for any remaining terms, their integral decrease rapidly as a function of v and the 
contribution of these terms approach zero in the limit R —>■ oo. Analogous argument 
can be repeated for the last integral in the rhs of Eq. (IB.6j) . only with the replacement 
a — > R/rj. Therefore, we have shown that the left hand side of Eq. ( IB. 61) approaches 
zero. Finally, the result of this analysis is that the only contribution coming from the 
annular region is the integral term in the rhs of Eq. (IB.lj) . Therefore, the summation in 
the Eq. ([3]) can be replaced by its integral. 
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